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Abstract. We study the perturbative power-series expansions of the eigen- 
values and eigenvectors of a general tridiagonal (Jacobi) matrix of dimension 
d. The (small) expansion parameters are the entries of the two diagonals of 
length d — 1 sandwiching the principal diagonal which gives the unperturbed 
spectrum. 

The solution is found explicitly in terms of multivariable (Horn-type) hy- 
pergeometric series in 3d — 5 variables in the generic case. To derive the result, 
we first rewrite the spectral problem for the Jacobi matrix as an equivalent 
system of algebraic equations which are then solved by the application of the 
multivariable Lagrange inversion formula. The corresponding Jacobi determi- 
nant is calculated explicitly. Explicit formulae are also found for any monomial 
composed of eigenvector's components. 
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1. Introduction 

The problem of solving algebraic equations by series expansions has a long history. 
In case of a single polynomial equation a major cornerstone is Birkeland's paper 
where a solution is given through an application of the Lagrange inversion formula 
|2| . Later on, an alternative approach was suggested by Mayr in who derived 
the relevant hypergeometric series as solutions of some PDE's satisfied by the zeros 
as functions of the coefficients of the polynomial. A modern interpretation in terms 
of GePfand-Kapranov-Zelevinsky hypergeometric functions E] can be found 
in jjj. The latter approach can also be applied to general systems of algebraic 
equations. 

The main goal of this paper is to derive a complete power-series solution of the 
spectral problem for a finite (d x d) Jacobi matrix M. We consider the off-diagonal 
matrix elements to be small, so that the whole problem looks like a perturbation 
of a diagonal matrix. After fixing a normalization of the eigenvector the problem 
is reduced to solving a system of d quadratic equations for d unknowns. Then we 
transform it into an equivalent larger system of special (Lagrange-form) 3d— 5 cubic 
equations, which are then inverted by the application of the multivariable Lagrange 
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inversion formula. The expansion of an arbitrary monomial of the components of the 
eigenvector is given explicitly in terms of multivariable (Horn-type) hypergeometric 
series of 3d — 5 variables. In the special case of an eigenvalue growing from a corner 
matrix element the number of expansion variables drops to 2d — 3. 
Consider a tridiagonal (Jacobi) matrix of order d 

71 a 2 



M = 



ih 



7k- 1 



o.k (3k 



(1.1) 



\ 7d-l a dJ 

and the corresponding eigenproblem MV — AV for the eigenvector V and the 
eigenvalue A. 

Assume that the off-diagonal elements (3k, Ik are the small parameters of the 
power expansion and that a.k are distinct. In the zeroth approximation (3k — 7fc = 
the matrix M is diagonal, its eigenvalues and eigenvectors being a& and V^ k \ 
respectively, where the components Vj of can be chosen as Vj® — Sjk- By 
continuity argument, for the small values of (3k-, 7& the eigenvalues A&, k = 1, . . . , d, 
are distinct and can be numbered in such a way that 

Afe = ctk + higher order terms. (1-2) 

We also choose to normalise the eigenvector V^ k ' 

MV (k) = Ak y(k) ( L3 ) 

by the condition 

Vt k) = 1 (1-4) 
for its k th component. Therefore, the remaining components must vanish in the 
zeroth approximation: 

Vj = + higher order terms, j ^= k. (1.5) 

The eigenvalue problem l|1.3(l together with the normalization condition (|1.4(l 
produces a system of d algebraic (quadratic) equations for the eigenvalue A/, and 
the free components vj k \ j ^ k of the eigenvector defining them as algebraic 
functions of the parameters a, (3, 7. The conditions i|1.2fl and 1|1.5|) fix uniquely the 
branches of the multivalued algebraic functions for small values of (3 and 7. 

The problem we solve in the present paper is to find an effective way to con- 
struct explicit expressions for the coefficients of the power series expansions for 

( k) 

the eigenvalues Afc and the components of the eigenvectors Vj '. The available 
literature on solving systems of algebraic equations by multivariate hypergeometric 
series (see a review in [7]) focuses mainly on solving the generic algebraic system. 
We are not aware of any detailed analysis of the particular systems arising from 
the Jacobi matrix spectral problem. The importance of the latter problem for nu- 
merous applications in mathematical physics and, in particular, in the theory of 
quantum integrable systems has been the main motivation of our study. Rather 
then using more modern approaches, in the present paper we follow the original 
idea of Birkeland and use a variant of the Lagrange inversion formula. 

We shall use the following variant of the multivariable Lagrange inversion the- 
orem jH| EI- Let boldface letters denote vectors £ = • . . ,£d), multi-indices 
q = (qi, . . . , qn) € Z D , and monomials £ q = . . . £|f . The inequality q > p is 
understood component- wise: qi > pi, Vi. Let [^ (J ]ft.($) denote the coefficient at £ q 
in the power series h(£). 
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Theorem 1.1. Let £ = (£ x , . . . ,£ D ), r\ = (771, . . . ,r) D ). Let ipi(r)), i = l,...,D, be 
formal power series in r\ such that tpi(0) =/= OVi. Then the system of D equations 

6 = -77, i = l,...,A (1.6) 

defines uniquely T]i(£), i = 1, . . . , D, as formal power series in £. 

In addition, let xiv) be a multiple Laurent series, that is x{ r i) r l p * s a power 
series for some p > 0. Then the Laurent series expansion for x{v(£)) 

x(i(0)= E ( L7 ) 

q>-p 

is given by the formula 

x q = \e]xHi)) = lv q Mv)v q (v) J(v), (1-8) 



where J is the Jacobian 



j =( ,,,[1_^| =< U«. 
a log J7 



(1.9) 



The analytic version of the Lagrange theorem |S] guarantees that if 1/^(77) and 
xi r l) r l p are analytic at then the Laurent series l|1.7l) converges in a punctured 
neighbourhood around £ = 0. 

When applying the Lagrange formula (|1.8J) . the major complication comes from 
the Jacobian J, which may be difficult to compute. Fortunately, for our particular 
problem the Jacobian can be calculated explicitly, in a relatively compact form. 

The paper is organised as follows. In section[2]we write down the set of quadratic 
equations defining the eigenvalue A and the components Vj of the eigenvector V, 
transform them into the form which is convenient for study, identify the combina- 
tions of the small parameters which serve as the expansion variables, and rewrite 
the equations again in the form which allows us to apply Lagrange's inversion for- 
mula. In section |3| we calculate an important component of Lagrange's formula: 
the Jacobian J. 

In section 0] we put together all the ingredients of the Lagrange formula and 
produce explicit expressions for all A^'s and vj^'s as finite sums of power series. 
The number of terms in the sum equals to the number of terms in the Jacobian J. 
All the power series are particular cases of a single universal Horn-type multivariable 
hypergeometric series which we denote $. In generic situation, the series $ depends 
on 3d — 5 expansion variables and d— 1 integer parameters. In section[S]we describe 
the simplification of our results for the special case k — 1 (or k = d) when the 
eigenvalue A stems from a corner of the matrix M . In section we examine a few 
low-dimensional examples illustrating the general results. The last section contains 
a discussion of possible applications and extensions of our result. 



2. Lagrange equations 

From now on we shall concentrate on studying a single eigenvalue A^ and the 
corresponding eigenvector V^ k \ for some fixed value of the index k (|1.2(l . We shall 
change our notation accordingly, to simplify the calculations. Set r = d — k and 
r = k — 1, so that 

d = r + r + l. (2.1) 

Let Afc = a/,. -I- A, Vj — Vk-j and ctj = ak-j for j = 1, . . . , k, and Vy = Vj-k 



and ctj = for j = k,...,d, so that V k — vq — vq and at — ao = a o- 

Respectively, let (3j = b^-j-i and a 3 = ak-j-i for j = 1, . . . , k — 1, and (3j — fij-k 
and ctj — ctj-k for j — k, . . . , d — 1. Without loss of generality we can set otk — 
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ciq = ao = 0. The normalisation condition (|1 .4|> implies that Vu = vq = Vq = 1. 
Since aj are assumed distinct, we have ai ^ and on ^ for i 0. 
As a result, the eigenvalue problem (|1 . 3f> takes the following form: 









\ 


/Vr\ 








h-i 






Vi 




Vi 




ci ai 






Vi 




Vl 




CO 


b 




1 


= A 


1 






c ai &i 




Vl 




Vl 






Ci-l 


a t bi 






Vi 


V 








\v r ) 




\v r J 



(2-2) 

The numbers r, r = 0, 1, . . . measure the distances of the selected diagonal element 
from the corners of the matrix. The cases r, r = 0, 1 are slightly special, we shall 
comment on them in due course, all other cases with r , r > 2 are generic. 
Expressing A from the 'central' ('zeroth') row as 



A = cqVx + b vi 



(2.3) 



and substituting it into the remaining rows we get the set of r + r = d — 1, see l]2.ip . 
quadratic equations 

Cj_iUi_i + aiVi + biV i+ i = c viVz + b viVi, i = l,...,r, (2.4a) 
CiV i+ i + aiVi + 6,_iUi_i = CQViVi + b viVi i = l,...,r, (2.4b) 

where we assume that vq = vq = 1, b r = 0, c% = 0. Formula l|2.3|) eliminates the 
eigenvalue A by expressing it in terms of the two eigenvector's components: vi and 
Vi. From now on, the variables Vi and Vi will be our only d — 1 unknowns. 

Note that equations l|2. 4J1 are invariant with respect to the 'involution' (rotation 
by 180 degrees around the central element) 



r <-> r, Vi <-> Vi, ai 



at, 



<-> Ci, 



(2.5) 



which we refer to as the tilde-symmetry. 

The next step is to rescale the variables Vi and Vi in order to identify the con- 
venient combinations of the expansion parameters. Setting bo = • • . = &r-i = 
and Co = ... = Cf-i = and denoting the corresponding values of Vi and as, 
respectively, vf and Sj 1 we get the equations 



Ci-iv^ + a^ = 



— — 7 ~-Q 







Ci-l o 
V 



i-ll 



(2.6a) 
(2.6b) 



which can be solved recursively yielding 

,.o _ (_\y c ° • • • Ci ~ 1 
ai . . . a t 

bo... bi-i 



~o 

V7 



i-iy 
i-iy 



ai 



(2.7a) 
(2.7b) 
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Let us rescale Vi and Vi by the formulae 



v % = v°iii = (-1)' C ° ' ■ ■ C% 1 m, i = l,...,r, (2.8a) 
a\...di 

v i =tfu i = (-l) i b °J-- bi - 1 u i , i = l,...,r (2.8b) 
a\ . . . di 



and rewrite equations 12.411 in terms of itj, 



fr c fr c „ frjQ . . . 

Ui = itiiti-- — itiUj H u i+1 , i = l,...,r, (2.9a) 

froco froc ~ bjCj „ . „ 

Ui = Uj_i - ^r- uiitj ^MiMi + ^c ttj+i, z = l,...,r, (2.9b) 

ami aiai ajdj+i 

where we assume uo = ^fl = 1, fr r = c r = 0. 

Equations l|2.9|l contain the parameters a, b, c, a, b, c in specific combinations, 
which are convenient to use as the expansion parameters. Note that bi and Cj enter 
only through the product biCi (same for b{ci). Introduce the variables Xi, yi, Zi 

6o co 



b c 

— ; 

fri+l C i+l 



i = 0,...,r-l, (2.10a) 

i = 0,...,r-l, (2.10b) 
i = 0,...,r-2, (2.10c) 



(altogether r + r + (r — 1) = 3r — 1 variables) as well as their tilde-analogs 

Xi = ^° c ° , i = 0, . . . ,r — 1, (2.11a) 

fro Co 



frj+lCi+1 



i = 0,...,f-l, (2.11b) 

i = 0,...,f-2, (2.11c) 

(altogether r + r + (r — 1) = 3r — 1 variables). The total number of variables is 
then D = 3r + 3F - 2 = 3d - 5. 

Equations l|2.9[) are simplified now to the form 

Ui = Ui-i - Xi-iuim - yi-xUiUi + i=l, ...,r, (2.12a) 

^ = - Ej-iSiM; - yi-xu\v,i + Zi-iu i+1 , i=l,...,r, (2.12b) 

where we assume uq = uq — 1, z r _i = z?-i = 0. Since fri,fri, Cj,Cj are small 
parameters, so are Xj, t/j, Zi,Xi,yi,Zi. In the zeroth approximation we have a set of 
binomial equations w° = w° = u^_ 1 . The iteration of equations l|2.12|l with 

the initial values u° = w° = 1 produces formal power series expansions of and 
in the variables (as, y, z, x, y, z) = 

The tilde-symmetry l|2.5|) for equations l|2.12l) takes the form 

r <-> r , Ut^Ut, Xi<r+Xi, yt^yi, zi^zi. ( 2 -13) 

Important remark: the variables £ are, generally speaking, not independent. 
Indeed, the 3d — 5 variables are bound by the relations 

^i = |^, * = 0,...,r-l, , = 0,..,F-1, (2.14) 
Vi x j fro co a i 

altogether r + F— 1 = d—2 independent relations, which leaves only 2d — 3 indepen- 
dent variables. Only in the special case r = (or r = 0), see end of this section, the 
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variables £ become independent. Nevertheless, when solving equations l|2.12() by 
power series, it is convenient to treat £ as the set of 3c? — 5 independent variables. 
The expressions (I2.10|) - (|2.11() for £ in terms of (a, b, c, a, b, c) can be substituted 
then into the resulting expansions at the very final stage. 

In order to find explicitly the coefficients of the power series for Ui and Uj we shall 
use the Lagrange inversion formula. In the notation of Theorem ll.il the expansion 
variables vector £ is composed of the 6 sets of variables £ = (x, y, z, x, y, z). To 
meet the theorem's premises we have to introduce also the vector r\ composed, 
respectively, of the 6 matching sets of expandable quantities 7] — (s, t, w, s, t, w): 



Si ~ Xi 


i = 0,...,r- 1; 


S 3 r 


° X j 


J = o, . . 


T — 1 


(2.15a) 


U ~ yi 


i = 0,. ..,r- 1; 




- Vj 


J = o, . . 


.,f-l, 


(2.15b) 


Wi <~ Zi 


i = 0,...,r-2; 


Wi r 




i = 0,.. 


.,f-2. 


(2.15c) 



Note that the variables (u, u) cannot be used directly as tj because they have 
nonzero limits for £ — > 0, therefore they have to be completed by small factors. 
Besides, the number of the variables (u, u) is only r + r = d — 1, so we need 2c? — 4 
extra variables to match the number D = 3d — 5 of expansion variables. Actually, it 
is more convenient to split the number D differently: D = 3c?— 5 = (2c?— 2) + (c?— 3), 
with the first set of 2c? — 2 variables for the larger system coming from only two 
variables, u\ and u\. Define Si, ti, s^, U as 

i = 0, ...,r-l, (2.16a) 
i = 0,...,f-l. (2.16b) 

The second set of d — 3 remaining variables w, w is defined as 

i = 0, ...,r-2, (2.17a) 

i = 0,...,r-2. (2.17b) 

It might seem easier to define Wi to be proportional to Ui rather than to the ratio 
of two u's with the step 2. However, a little experimenting shows that our choice 
of Wi leeds to simpler, factorised expressions for the functions <Pi(rj) in (|1.6(l . 

From Mo = uq = 1 it follows that 

w Q = z a u 2 , (2.18a) 
w = z u 2 . (2.18b) 

Solving equations (|2.17(l - (|2.18|l recursively we obtain the expressions for u 2 , . ■ ■ , u r 
in terms of Wj and for u 2 , ■ ■ ■ , in terms of Wj : 



Si = XiUi, U = yiUi, 

Si = XiUi, U = yiUi, 



Uj+2 
i ■> 

Ui 

Ui+2 
i ~ j 

u,. 



w 2j - 2 ...w W 2 j-1 ■ ■ ■ WlSo 

U 2 j — , u 2 j+i — i 

Z 2 j-2 ■ ■ ■ Z Z 2 j-1 ■ ■ ■ ZiXq 



(2.19a) 
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Proposition 2.1. The set of equations \2.12\) . $2.16)) and \2.11\j is equivalent to 
the set of Lagrange-type 11. b)) equations 

Si t{ 

' Hi 1 " ") • • • ) ' A 7 j 



x% = —, Vi = —, i = 0, ...,r- 1, zj = — , j = 0, ...,r-2, 



(2.20a) 



Si — tl ■ r\ — -1 — j 



Xi = — , Vi = —i i = Q, ...,r-l, zj = sr- , j = 0, ...,r-2, 

(2.20b) 



Jo = • • ■ = /r-i = ,g = • • ■ = = — — , (2.21a) 

1 + s + t 



hi = + v -o,..,,-2 k^.o), (2.21b) 

(1 + Si + ti){\ + s i+ i + t t+ i) 

7 7 1 + ™° m oo ^ 

,/o = • • ■ = J?-i = .go = • • ■ = ffr-i = — , (2.22a) 

1 + so + to 

'■ - + i = ,... ) r-2 (£^0). (2.22b) 



{l + Si + ti)(l + S i+1 +ti + i) 

In the notation of Theorem ll.il we have 

£ = (x,y,z,x,y,z), r) = (s, t, w, s, t, w), <p — (/, g, h, f,g,h). (2.23) 

Proof. Consider equation (|2. 12a|) for Uj+i: 

u i+ i =m — XiUxm+i - yiUiU i+ x + ZiU i+2 . (2.24) 
Note that ZiUi+2 — WiUi by virtue of (|2.17a|l . Having rearranged the terms we get 

+ XiUi + t/jtti) = Mi(l + Wi). (2.25) 
From (|2.16a|l it follows that XiU\ = Sj and y^ui = i^. Therefore, 

JiL = l±it±*t. (2.26) 

ltj+1 1 + Wi 

Multiplying the equalities (|2.26() for the index i and i + 1 we get 

Uj _ (1 + Sj + tj)(l + s t+ i +tj+i) 
Ui+2 (1 + Wj)(l + Wi+l) 

It remains to replace Ui/ui+2 with Zi/wi from the equality (|2.17a|) . and we obtain 
the equality of the form z% — Wi/hi, see (|2.2()af) . where 

(1 + ^(1+^) i = ,... I r-2 > (2.28) 

{l + s ^ +t l )(l + s l+1 +t t+1 ) , y J 

and it is assumed that w T -\ = 0. 

The special case of equation (|2.12al) for i = 1, 

ui = 1 - x iti — youxux + z u 2 , (2.29) 

has to be treated separately. Similarly to the general case, we replace Z0U2 with wq 
using (|2.18a[l . and substitute u± = sq/xq from H2.16a|l and U\ — 'sq/xq from 12.16b|l . 
Then we replace yoSo with xoto, using both l|2. 16a|) and l|2.16b|l . As a result, the 
equation takes the form 

s (l + s Q + t ) = x Q (l + w ) (2.30) 



(2.27) 
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or, equivalently, the form x$ — S0//0, see (|2.20a|l for i = 0, where 

A = ^-. (2.31) 
1 + so + 

From (j2.16a.fl and l|2.16b|l it also follows that 

f i = S -± = u 1 = %=g i = S -±=f . (2.32) 

Xi yi x 

The remaining half of the equations is obtained by the tilde-symmetry (12.1311 . ■ 
As was said before, the cases of a corner eigenvalue, i.e. when r = or r = 0, 
and the next one, i.e. when r = 1 or r = 1, are slightly special. When one applies 
the formulae from the above theorem describing the generic case, i.e. when both 
r, r > 2, to such cases one has to remember that 

• for a corner eigenvalue, say r = 0, there are no tilded variables and there 
are also no variables y,t,g, so that in this case we only have 2d — 3 = 
2(r + 1) — 3 = 2r — 1 variables, therefore one must disregard (|2.22|l entirely 
and set tj = 0, j = 0, . . . , r - 1, in I2.21[) : 

• for the immediate next eigenvalue, say r = 1, all variables are present but 
w, so that one must set wo = in (|2.22a(l and disregard (|2.22b(l entirely, 
there are 3d — 5 = 3(r + 2) — 5 = 3r + 1 variables in this case. 

All these modifications are easily seen from the expressions (I2.10(l - (|2.11|l of the small 
parameters (x, y, z, x, y, z) in terms of the initial small parameters (b, c,b,c). For 
a more detailed study of the special cases see section GD 

3. Jacobian 

In this section only, we use the ordering of variables which is different from the one 
used above: £ = (x, y, x, y,z,z),r] — (s, t, s, t,w,w), tp = (/, g, f, g, h, h). 

Theorem 3.1. The Jacobian J defined by \1.9\) can be expressed as 

J = SS-TT, (3.1) 



where 



■l^l + s j+tj U Wk ' T -W^l + s j + t j 11 ^ (3 ' 2) 

j=0 J J k=0 I j=0 J J fc=0 



and 



r-2 



W=\[(l + w k ). (3.3) 

k=0 

Here we adopt the following agreement: whenever a sum has the upper limit 
smaller than the lower one then its value should be taken as 0. Also, any product 
in a similar case should be taken as 1. The tildes in l|3.1|l refer to replacing Sj, 
tj, Wj and r by their tilded versions. As always, we assume that Wf-i = and 
w r -i = (see Proposition ^. 1[1 . 

Note that the number of terms in J (13.1(1 is rr + (r + l)(r + 1) = 2rr + r + r + 1 = 
2dk - 2k 2 + 2k - d (recall that d = r + r + 1 and k = r + l). 

Proof. Consider the rows of the matrix J = S-jh — — ir^'- 

— — — — ^(3 4) 
1 — — 
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(o,o,ia...,o,...,o,...,o,...,o,^, t ^a...,o,. 



(o,... ) o,i J o,...,o,...,o,..., IT |^ )IT | T j r ,o,...,o,... 
(o,... ,0,0,1,0,. ..,o,...,o,...,o, TT ^, T ^,o,...,o,. 



( n a i -i a n n ^ n i p*} 

\ u > ■ ■ ■ ' l + so+to ' ■ ■ ■ ' ^ 1 + so+to ' 1+so+to ' ■ ■ ■ ' ' ■ ■ ■ ' ' ■ ■ • ' 1+so+to ' ' ' ' 7 ^ ' 

(o I). 1.1) 

(o,...,o,...,o, o,i,o,... ,o,...,o,...,o, TT ^, TT f T ^,o,. 



i n *s! *o i n n ^ n i 7) 

\ u > ■ ■ ■ ' l+so+io ' " ' ' l+so+to ' ■ ' ■ ' ' ' ■ ' ■ ' ' • ■ ' ' l+so+to ' u ' ■ ■ 7 I 

(0,...,0,...,0,...,0,l,0,...,0,..., TT i^ 7 , TT i Tir ,0,... 

(o,...,o,...,o,...,o, 0,1,0,... ,o,...,o, TT ^, TT ^,o,. 



(- T |^,...,0,...,0,...,- T ^,..., T ^,0,...,0,...) (3.8) 

(o,...,o ) ... > o 1 ...,o,... > - T ^ r , I ^ sr> o J ...,o J .. 

(o,. .. ,0,. .. ,0,. .. ,0,. .. ,0,-^,^,0,. ..,0,. 



(o,.. . ,0,. . .,0,.. . ,0,. . .,0,. .. ,-^-,^-,0,. . 

(o,...,o,...,o,...,o,...,o,... ,0,-^,^,0,. 



First lines of each of the four first sub-sets above, which correspond to the variables 
So EH , t ((3.5(1 . so (13.6(1 and t ((3.7(1 . can be replaced by simpler versions (without 
changing the determinant) by adding a multiple of the row associated with wq ((3.8(1 
or wq ((3.9|) . They become as follows: 

(l s (i+w ) 1 

(o 10 o o * o(1+ y o) 1 
(o n i (i n n so ( l+^o ) n A 

^U, . . . , U, . . . , 1, u, . . . , u, . . . , u, . . . , ( 1+So+to )„ , u, . . .j 
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t a (l+w ) 



(o,...,o,...,o,...,i,o,...,o,..., 



(l+s +t )w 



,0,...) 



After the above replacement, the (3r + 3r — 2) x (3r + 3r — 2) matrix J' 
— ^- acquires the following block-matrix form: 



J = ^(2?+2r)x(2? + 2r) ^ ^ J _ ^ ^ = _ ^ 



where the matrices B, C, D are 



/ S (l + M>o) 

' (l+3q+to)i"o ' 
si 

l+3i+ti ' 



0,..., 

to(l+w ) 
(i + so+t )w ' 
ti 

1 + si+tl ' 

o, 



0,.... 

0,.... 
0,.... 
0,... : 

0,.... 
0,.... 
0,.... 
0,.... 

V o,.... 



si 



1+sj+ti 

S2 

l+5 2 +*2 ' 



o, 

tl 



1+sj+ti 

l + S2+t 2 ' 



o, 

o, 
o, 
o, 

o, 
o, 
o. 

o, 



l + S2+*2 ' 
1 + Sf-l+tf - 

o, 
o, 

*2 

l + S 2 +*2 ' 



l+s ? _i+*?- 

o, 
o, 
o, 

o, 
o, 
o, 
o, 



0,..., 
0,..., 
0,..., 

0,..., 
0,..., 
0,..., 
0,..., 

0,..., 

S (l + W ) 

(l + s„+t )w ' 



1 + si+tl ' 

o, 



0,..., 

t (l+ui ) 

(l + s +t )tuo ' 
tl 

1 + si+tl ' 

o, 



o, 
o, 
o, 

o, 
o, 
o, 
o, 

o, 
o, 

«1 



1+si+ti ■ 
si 

l + S 2 +«2 • 

















S2 



l + S 2 +*2 



o, 


Sr-1 


l + ST ._l+i,._l 







tl 





1 + si+tl ' 


*2 


tl 


l + S2+*2 ' 


l + S 2 +*2 



o. 







l + Sr—l+ir 



l+w ' ■ ■ ■ ' 

0,..., 



0,..., 
0,..., 



0,... : 
0,.... 



m_ \ 

l+w a ' • • • \ 

0,... 



c = 



0,..., 
0,..., 
0,..., 



0,..., 

WQ 



1 + WQ ' • • 

0,..., 



0,.. 



? " " " 5 

WQ 



1 + WQ ' • 

0,... : 



0,.. 
0, .. 
0,.. 



V o,..., 



0,..., 



0,..., 0,... / 
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1+wo ' 

1+Wl 



0, 

1 

I+tui 



0,.. 
0,.. 



0. 
0. 



0,.. 
0,.. 



o \ 




D 



0,... 


fi>r-2 


1 


l+S>r-2 ' 




0,... 


o, 


o, 


0,... 


o, 


0, 



0. 
1 

1+wo ' 

_ Wi 



0,... 

0, 

1 

1+11)1 



(3.10) 



0,. 



0, 



Wr-2 



1+tt), — 2 ' 1+U) r _ 2 / 



Therefore, the calculation of the Jacobian has been reduced to finding an explicit 
formula for the determinant of a smaller, (r + r — 2) x (r + r — 2), matrix D — CB. 
Subtraction of the product CB does only change two rows of the triangular matrix 
D, those associated with the variables wo and wq- They now become 

tig ( 1 , S (l + S ) . Si 

l+w \w (l+s +t )w l+st+ti' 

Si j S 2 ^ ^r-2_ I 

l+3l+tl l+S 2 + t 2 ' " ' ' l+S ? _ 2 +t?-2 l+Sy_l+t^_l ' 

to(l+mo) . tx 
(l+s +to)w ^ 1+si+ti ' 

tl I t 2 tr-2 I *r— 1 A 

1+Sl+tl T " l + S 2 +t 2 ' ' " ' l+S r _2+tr-2 l+S r _l+t r _l ^ 



and 



"0 [ *o( 

1+UJQ \(1 + S< 



1+yo) 

So+to)wo 

h i 



1+Sl+tl 
*2 



1+Sl+tl l + s 2 +t 2 ' 
1 , s (l+tu ) _^ si 
Wo (l + s +to)w 

si 



l + Sf.-2+tf 



1+si+ti 



1 + si+ti • 

S2 

l + S 2 +«2 : 



l+S r _ 2 +tr 



l + S ? _ 1 +t f ; 



l + Sr ._l+t r _l 



respectively. The rest of the rows in the matrix D—CB are the same as in the matrix 
D above, that is they all have only two non-zero elements: one is on the diagonal 
and the other one is to the left of it. By adding multiples of the columns, without 
changing the determinant, we can arrange that all those non-diagonal entries in the 
selected (r — 2) + (r — 2) rows vanish, allowing us to compute the determinant by 
reducing it down to a 2 x 2 determinant. 

Indeed, multiply the last column of the matrix D — CB (look at formula l|3.10(l ) 
by w r -2 and add the result to the penultimate column. The last row has now got 
only one (diagonal) element. Hence, we keep the factor — - and reduce the 
determinant to its minor, removing the last column and the last row. By repeating 
this process, we shall end up with the expression 



wo wo (SS — TT) 
(1 + w ) •••(! + u) ? _ 2 )(l + w ) •••(! + iu r - 2) ' 
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where 

r— 2 

o _ 1 , s (l+m ) , ai(l+mi) , Sjtl+tuj)!]^ 1 ! ">fc , gr-1 UlZl w k 

° w (l + s a +t )w„ ~T 1 + si+ti f L 1 + Sj+tj l + s r _i+t r _i ' 

3=2 

r-2 _ ._ 

§ = _L + So(i+a ) , si(i+ai) , sj(i+«jj)nfc = l ™fc , s r-i ]lfc = i w k 

3=2 

r-2 

rf t (l+W ) ti(l + mi) V^ t J (l+<gj)I]U'""= I tr-lUlJ^k 

(l+s +t )w + "I" 1+Sj+ij "f" l+Sr-l+tr-l ' 

r-2 _ 

(l+*o+*o)'5)o 1 + si+ti A-^t 1+sj+tj l+sp_!+tp_i ' 

J=2 

This is apparently equivalent to the statement of the Theorem. ■ 



4. Hypergeometric series 

Now we have all the ingredients for the right-hand side of the Lagrange formula i 
and can calculate the expansion coefficients. In this section we return to the initial 
ordering £ = (x, y, z, x, y, z), r] = (s,t,w,s,t,w) and <p = (f,g,h,f,g,h). 

Let us choose the function xiv) m I 1 -ty to be a Laurent monomial x(v) = rf = 
s m ' t n ' wP's™'!™' w& , q' eZ D = Z 3d - 5 . 

Define the step function <xy as 

?' ' ' ■ (4-1) 

1, i>3 y ' 

We shall use the binomial, trinomial and quadrinomial coefficients defined for inte- 
ger m, n,p as 

a )^{x m ](l + xT, ( a ) = [ x ™y»](l + x + y)*, (4.2a) 
ml \m,nj 



x m y n z p ](l + x + y + z) a . (4.2b) 

ats are evaluated as 
a\ a(a — 1) . . . (a — m + 1) / a \ a(a — 1) . . . (a — m — n + 1) 



Note that the multinomial coefficients are evaluated as 



ml ml \m,n 



a \ a(a — 1) . . . (a — m — n — p + 1) 



y m,n,pJ m!n!p! 
for m > 0, n > 0, and p > and vanish if m < 0, or n < 0, or p < 0. 
Let |m| = mo + . . . + m r _i etc., and set 



(4.3a) 
(4.3b) 



p_ x = |m| + , p_ x = |m| + \n\ . (4.4) 

We shall also assume that p r -x = Pf-i = 0. In the rest of the paper we shall 
frequently present formulae for the untilded quantities only, assuming, unless oth- 
erwise stated, that the tilded versions are obtained by tilde-symmetry. 

Theorem 4.1. The expansion of the monomial r] q in £ q is given by 

q>q' 
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where 

U q . q = [€ q }ri q ' = [x m y n z p x^y h zP]s m 't n 'w p '^'t Tl 'wP' 

r — 1 r—1 r— 1 r — 1 

= E E%-EE% (4.6) 



i=— lt=— 1 i=0 7=0 



J=o V'"./ ~ '"' - 5 ij , ".; "' / V /'.. - />' - o-y / 

gi ; = TtV 1 _ ~ s *i ^ ( Pi-i + /'• + s *i 1 \ , . _, 7 , , , 

J ' 1 '" 7=o V m J ~ m i' n i ~ n 'i ~ "W V Pj -Pj-Vij * 



(and respective tilded versions). 



Proof. Using the definitions 1(2.21(1 and 12.22(1 we get for q — (m, n,p, m, n,p) : 



<p*(r,) = f m g n hPf m g n h p 

— JO •••J) — 1 tfo •••»r-l "0 ■■■"■r-2 

x /0 3o ■■■9r-i n ■■■ n r-2 



n t 



1 + s + t J \ 1 *- V 1 + s, ■ + tj 

, - , ,„i|+|n|+p /r-X / , . ~ \W-i+W N 



H r 



1 + so + to/ V7=i V 1 + S J 

where we assume w r _i = iuy-i = £> r -x = i?f-i = 0. 

Using the shorthand notation (|4.4(l we can write down <p q (ri) in the compact 
form 

, vpm+pA /V-i / -, . ~ 

1 + W,- \ TT / 1 + W 7 - \ 



Using the step function <7y 1(4. 1(1 we rewrite the expressions 1(3.2(1 for the ingre- 
dients of the Jacobian J in the following equivalent form: 



r -E!!7r 



W 3 V 



5 = V TT j J 

<~i|=o (! + w j) l ~ Sii (! + fl J + *i) aw ' 

Substituting the above expressions for y> 9 and J into the right-hand side of 
the Lagrange formula p. 8(1 we get 

r—1 r—1 r-lf-1 

J -EE ^ -EE gi &> ( 4 - 8 ) 
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where 



r-l 

X 

3=0 
r-l 



ft v l + sj + tj) v i + sj + tj j (i + Wj y- 5 *i 

x(l + Sj +tj) - p ' -4 « (1 + Wj +w +5 « - 1 , (4.9) 



^ = II af'tf^V^^ (1 + 8j + tj )- p i-^- 5 ^ (1 + Wj )P3-i+Pi+^-l (4.10) 

and J-" 1 , Q L are their tilded versions. In the above formulae it is assumed that 

<ri 1,r -' = o° = i. 

It remains to take the coefficient of the resulting expression for xy 9 J at the 
monomial r\ q . Using the binomial and trinomial expansions (|4.2al) we get the 
expressions (|4.7a|) and (|4.7b|) for the expansion coefficients F l q ' q = [s m t n w p ]T l and 

Q l , = [s m t n w p ]Q l . The final expression ({4.6(1 for [£ q ]p q follows then immediately. 
■ 

As formula (|4.(il) shows, the monomial r\ q expands into £ q as a finite sum (the 
number of the terms equals to the number of the terms in the Jacobian J) of 
Laurent series. These Laurent series have the uniform structure and can be in fact 
expressed in terms of a single standard series. 

Let us introduce the function $(£; /i., ji) of D = 'id— 5 complex parameters £ and 
depending on two integer vectors /i. = (yito, . . . , // r -i) G 27 and Jl — (juo, . . . , /ly-i) € 
Z r , altogether r + r = d — 1 integer parameters: 

n\tq 



3(fc,i,£) = £(-l)M+M+lftl+l«l € 

q>0 

T— 1 / \ f— 1 



X 



■Q f Pi + m 3 + n 3 + Pj-1 + V] - 1 \ -Q 



in j j 5 



+ m j + n i + Pj'-i + — 1 



(4.11) 



Note that, despite the appearance, <&(£; fi, j5) does not factorise into tilde and 
no-tilde factors because of p-% = \m\ + \n\ 1)4.4(1 . 

For (Ltj > 1 the coefficients of the series ((4.11(1 can be written down in terms 
of the Pochhammer symbols (a) m = a(a + 1) . . . (a + m — 1): 

m!n!p!Tn!n!p! 

q>0 J ' 

x J~J (A'jQmj+Kj+Pj +Pj-i J-J j ) fhj +nj +pj +pj _ i ^ ^ 

where m! = mo! . . . m r _i! etc. 

Formula 1(4.12(1 characterizes $ as a generalised Appel-Horn-type series: the 
ratios of adjacent coefficients of the power series are rational functions of the indices 
|10|. Note that $(0) = 1, and the series converges in a neighbourhood of |1(J|. 
The function $ can also be viewed as an „4-hypergeometric function in the sense of 
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Gelfand-Kapranov-Zelevinsky |SJ |B] — we plan to elaborate on this remark in a 
subsequent paper. 

Theorem 4.2. The expansion of the monomial r\ q in £ q can be expressed in terms 
of the function $(£; /x, ju) as follows: 



rf' = £i' ^ t' ' «H{: 7T + v\* + u X ) 

Vi=-1T=-1 
r— 1 r— 1 \ 

-EE^^^+^'^+^L ( 4 - 13 ) 

i=0 1=0 / 

i = -1 

€"= | | x^z^ = { x , i = (4.14) 

i > 




whe 



i = -1 

C" IK :J { 2/0, * = (4.15) 

i > 

7r J = P ;_ 1 +^, y) = 2<T iij - U j = 0,...,r-l (4.16) 
and respectively for the tilded versions. In \4-lb\) it is assumed thatp'_ 1 = \m'\ + \n'\. 

Proof. Note that the sum over q > q' in 14.5JI can be safely replaced with the 
sum over q £ Z D , since the multinomial coefficients in (|4. 7|) vanish for the negative 
values of the lower indices and thus automatically select the correct limits of sum- 
mation. Then from H4.5(l and (|4.6J) we can express r\ q as the finite sum, each term 
corresponding to a term in the expansion (|3.1|l of the Jacobian J: 

r— 1 r— 1 r-lr-l 



/=EE ^)-EE G ^o ( 4 - 17 ) 

i=-l t=-l i=0 Z=Q 

= E ^V^*> Gj(0 = £ fifr,^,**. (4-18) 

Consider the single term F^, (£) of the above sum. From 14.7a(l it follows that the 
summation can be restricted to the values of the indices for which the multinomial 
coefficients do not vanish: nij > m'j + Sij, rij > rip pj > Pj + fXy, resp. mnp. Let us 
shift the summation indices: rrij := rrij + rn'j + rij := rij +n'j, pj := pj +p' rj +o~ij, 
etc. so that the summation runs over q > 0: 

if (€) = £«V <+h J2 P'/k^' (4-19) 

q>0 

r = TT f-Pi-i ~ Pj ~ P'j-i - Pj - 2(T ^-i 
q ' q 11 \ rrij, rij 



3=0 

r Pi-l + Pj + P'j-i +Pj + 2( ?i,j 
Pj 

In doing so we have used formulae 14.7a|) . (|4.14|) . and the identity Sij + a v 

-l 



(4.20) 



(Ji.j-i ■ Note that shifting rrij in p_i = \m\ + |n| produces the term X^=o = a i,-i 
which is consistent with the convention p'_ x = \m'\ + |n'|. 
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Applying the identity 



-aWa-l\ fa + m + n-V 

to, nj \ p J \ m, n, p 



to l|4.20[l we get another expression for ¥ l , 



r-1 



p f = ^\m\+H TT ( m J + n 3 + Pd-i + Pi + Pj-i + P'j + 2 °U-i ~ A (4 22) 

Repeating the same steps for FS and G^/(£) and collecting together the ob- 
tained expressions we identify the result as H4.13fl ■ 

The expansions we are ultimately interested in are those of the monomials in the 
components Vj and Vj of the eigenvector V <|2.2|) . or of their rescaled versions: Uj 
and Uj (|2.8|) . 

Proposition 4.1. Let u k = (u kl , . . . , u^), keZ r , and resp. for u k . The expan- 
sion of the monomial u k u k in £ q is given by the formula 

r— 1 r—1 _ 

u k u k = £ fl+h $(£;7r + //,ff + 

i=-lT=-l 
r—1 r—1 



where £ 0i , f are same as in Theorem \4- 6 J\ and 

■Kj = k j+ i + . . . + k r , j = 0, . . . , r - 1, (4.24) 

Proof. Using formulae l|2.19|l we can express Uj as monomials in x$, soj z and w - 
Thus 

u k = x -m' z -p' s m> w p> ( 4 

where 

mo = fci + fc 3 + . . . , m[ = ... = m! r _ x = 0, (4.26) 

P'lj = k 2j+2 + k 2 j+4 + ■■■, P 2j+1 = k 2 j+3 + k 2 j+5 + ■■■ (4.27) 

(similarly for the tilded variables) . Substituting these values of p'j into the formula 
ttj = p' j _ 1 + p'j <|4.16[) we get (I4.24|) . Note that p'_ x = \m'\ = m' . It remains to 
apply Theorem l4.2l to the monomial (|4.25l) . ■ 
To obtain the expansion of the monomial v k v k in bebe one needs to express the 
components Vi of the eigenvector V in terms of uu using (|2.8|l . then apply Proposi- 
tion ^3 an( l finally using (|2.1()|l and 1|2.11|) . to express £ in terms of bebc. To get 
the expansions involving the eigenvalue A one has to use (12. In particular, the 
resulting expansions for ui or ui give an explicit solution for iterations of equations 

In this paper we study only formal power series expansions. The analytic version 
of the Lagrange theorem guarantees that our power series converge in an open 
neighbourhood around £ = 0. The precise description of the convergence domain 
of a multivariate power series is, however, a difficult task, see for example |6l I1U|. 
and we leave it to further study. 
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5. Corner eigenvalue 



As mentioned in the end of section the cases of the corner eigenvalue r = 
(r = 0) have some peculiarities. Let us consider the case r = in more detail (the 
case r = can be treated by the tilde-symmetry). We have 

/0 b \ 

c a-y 6i 



V 







/1\ 


Vl 




Vi 




= A 




Vi 




Vi 


\v r J 




\Vr) 



(5.1) 



The tilded variables are absent. Besides, the sequences yi, resp. ti are also absent 
The total set of variables £ = (x, z) has the cardinality (r)+(r— 1) = 2r— 1 = 2d— 3 
and all the expansion variables are independent, in contrast with the generic case 
Respectively, tj — (s, w) and ip — (/, h), where 

l + W (1 +Wi)(l + w i+1 ) 



fa 



fr-l 



i = 0, 



1 + V (l + Si )(l + Si+l) 

The expression for the jacobian J simplifies to 

J = s = fi+E n «*) 3d + ^r 1 



j=0 



fe=0 



fc=0 



and contains thus only r + 1 = d terms. The analog of formula (|4.8|l is 

r-l 



where 



rn'-+8ij p'+Cj. 



(l + s J -)- ft '- 1 ^'- 5w (l + ^) 



j=0 



q q' 



The expression 1)4.6( 1 for the expansion coefficient becomes 

r-l 

H q , q ee = [a; m zP] S ™'«,P' = £ T\ 

where 

r-l 

^„ = [s m ™ p ]r = n 



-Pj-i-pj 



,_„ x mj-m'j-Sij 



r-2 

n 

J=0 



r-2. 
(5.2) 

(5.3) 

(5.4) 
(5.5) 

(5.6) 
, (5.7) 



whith p_i = | ?7i | . 

The expansion of the monomial u k in x m z p is given by the analog of formula 

r-l 

ufe = E l fl *(^,2;7r + ^) (5.8) 

i=-l 

where £ , f l are same as in Theorem 14.21 7r is given by l)4.24[l . and 3>(x,z;fi) is 
defined as the series 



$(x,*;m)= E (-^)^ P I1 

m,p>0 J— 



+ TO i + Pj-i +Pj ~ 1 



(5.9) 
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or, for (Xj > 1, 



Hx,z;„)= X) ^n%^. (5.10) 
' m!p! J -- L (Ui)v , 

6. Examples 

In this section wc illustrate our general results with a few low-dimensional examples. 
6.1 The simplest case is d = 2 (r = 1, r = 0). The spectral problem is 

b " Ul ^ a^y (6.D 



s c aiy yuiy \yi y 

After eliminating the eigenvalue A = boVi and introducing the rescaled variables 
U = u\ = —viCi/cq, x = xq — boco/al we get for the single unknown variable u the 
single quadratic equation 

1-u-xu 2 = 0. (6.2) 

The branch we are studying is selected by the condition u — > 1 for x — > 0, and 
the corresponding solution is 



yi + 4a:-i 

U W = ^ ■ 6.3) 

2x 

In terms of the variable s = xu we get the equivalent Lagrange-type equation 

x = s(l + s) = L f = T^—- (6-4) 
/ 1 + s 

The corresponding Jacobian J has two terms: 

s df s , 

J = l--± = l + -—. 6.5 
/as 1 + s 

The expansion of the monomial £ Z given by Proposition ^. II is 

u fc = $(.x;fc) +a;$(a;;fc + 2), (6.6) 

where 



v-^ , > m ( li + 2to — 1\ 1 / vl + 4x — 1 \ 



m>0 

For /i > 1 the function $(a;; /i) can be expressed in terms of the Gauss hyperge- 
ometric series: 

As a matter of fact, for k > the two hypergeometric series in ()6.6|l can be 
summed into a single series 

6.2 Three-dimensional case, corner eigenvalue: d = 3 (r = 2, r = 0). The 
spectral problem 

b 0\ /1\ /l 

c oi hi Ui = A U | , ((i.lih 
ci a 2 / V v 2/ \«2> 
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after the substitutions A = boVi, v\ = — u\Co/a\, V2 = U2CQC\/a\a2 is reduced to 
the pair of quadratic equations for the variables u\, U2 

1 — Mi — x$u\ + zqU2 = 0, (6.11a) 

ui — ii2 — x\U\U2 — 0, (6.11b) 

where the expansion variables are 

boc b c feici 

x = — xi = , z = , (6-12) 

a| aia2 a\a,2 

and we choose the branch m, 112 — > 1 as xo, x\, zq — > 0. 

The equivalent Lagrange-type equations for the variables sq = xqUi, s% = X\U\, 
wq = Z0U2 are 

so 9i w 

x = —, Xi = —, z = —, (6.13) 
Jo Ji "0 

where 

1 + Wq 1 + wq 

Jo = fi = -j- , h = 77- rrr- r- (6.14) 

l + s (l + s )(l + si) 

The three-term Jacobian J equals 

J=^— + ^ + - ™ (6.15) 

1+Wq 1 + So (l + «i)(l + 1tfo) 

By formula (|5.8|l . the monomial u^u^ 2 expands in xq, Xi, zo as 

uj 1 ^ 2 = $(xo, xi,zo; &i + ^2, fe) 

+ x $(xo,xi,2o;fei +fe + 2,fc 2 ) 

+ Xiz $(x , xi, 2 ; fei + + 2, k 2 + 2), (6.16) 

where 



$(x ,x 1 ,^o;Mo,Mi)= E (-^o) mo (-xi) mi ^ 1 

mo,mi,po>0 

Vo + 2m + mi + po - 1\ ( fii + mi + p Q - 1 
m , po A mi 



or, for /i ,Mi > lj 
$(xo,xi,z ;mo,Mi) = X] 



(6.17) 

y/tl / 

(-X ) m ° (-Xi)™ 1 Zg" 1 (po)2m +m 1 +p (Pl)mi+p 



m „, mi , Po >o mo!mi!p " ! (Mo) mo+mi (Mi) P c 

(6.18) 

6.3 Three-dimensional case, middle eigenvalue: d= 3 (r = r = 1). The spectral 
problem 

cti 60 









.3 


-1 











c 6 I I 1 I = A I 1 I (6.19) 

C Oly 

after the substitutions A = Co«i + ^o^ii u i = —uco/ai, Vi = —ubo/ai is reduced to 
the pair of quadratic equations for the variables u, u 

1 — u — xv? — yuu = 0, (6.20a) 

1 — u — xu 2 — yuu = 0, (6.20b) 
where the expansion variables are 

boco &0C0 ~ b c „ & Co ,„ 01 \ 

» = — 2-, y=~ — , x = ^2-, y = ^-, (6.21) 

aiai a| aiai 

and the chosen branch is u, u — > 1 &s x,y,x,y —> 0. Note that the variables x,y,x,y 
are bound by the single relation irx = yj/. 
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The equivalent Lagrange-type equations for the variables s = xu, s = xu, t = yu, 
t = yu are 

S t „ 5 _ t . . 

x = -, y = -, y=~, (6-22) 

f 9 f 9 



where 



f = 9 = T T— t, 1=9= - ; i ; ~ ■ (6.23) 
1 + s + 1 1 + s + < 



The five-term Jacobian J equals 

s * 

J = 1+ ^ 7 + 

l+s+t 1+s+t 

sis tt 



(6.24) 



(l+s + t)(l + s + t) (l + s + t)(l + s + t) 
By formula (|4.23(l . the monomial u k u k expands in £ = (x,y,x,y) as 

u k u k = $(£; fc, fc) + ; k + 2,k) + fc, fc + 2) (6.25) 

+ a:x$(£; k + 2, fc + 2) - fc + 2, fc + 2), (6.26) 

where, for q = (m, n, m, n), 



- £(-aO TO (-i/) n (-2r(-y) n 

q>0 

/ /i + 2m + n + n — 1\ //I + 2m + n + n — 1 
m, n J \ in, n 



(6.27) 



or, for /i, /i > 1 

^ m!n!m!n! (|u)m+K(£0m+n 



q>0 



7. Discussion 

In the present paper we have solved a problem of much physical interest. The 
spectra of finite Jacobi matrices appear in many applications: from orthogonal 
polynomials and nearest-neighbours interaction models to solvable models of quan- 
tum mechanics (Lame polynomials and Bcthe ansatz). 

What is left for further study are the questions of convergence domains, differen- 
tial equations and integral representations for the obtained hypergeometric series, 
as well as their relation to A- hypergeometric functions introduced in 4,56. The 
approach used in our work can be generalised to multiparameter spectral problems 

Explicit perturbative solutions have also appeared in a different context in the 
works by Edwin Langmann |12l I13| . notably for the multi-dimensional spectral 
problems related to the Calogero-Sutherland and elliptic Calogero-Moser systems. 
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